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ABSTRACT 

Context. We explore the physics behind one of the brightest radio afterglows ever, GRB 030329, at late times when the jet is non- 
relativistic. 

Aims. We determine the physical parameters of the blast wave and its surroundings, in particular the index of the electron energy 
distribution, the energy of the blast wave, and the density (structure) of the circumburst medium. We then compare our results with 
those from image size measurements. 

Methods. We observed the GRB 030329 radio afterglow with the Westerbork Synthesis Radio Telescope and the Giant Metrewave 
Radio Telescope at frequencies from 325 MHz to 8.4 GHz, spanning a time range of 268-1 128 days after the burst. We modeled all 
the available radio data and derived the physical parameters. 

Results. The index of the electron energy distribution is p = 2. 1, the circumburst medium is homogeneous, and the transition to the 
non-relativistic phase happens at f^R ~ 80 days. The energy of the blast wave and density of the surrounding medium are comparable 
to previous findings. 

Conclusions. Our findings indicate that the blast wave is roughly spherical at fNR, and they agree with the implications from the VLBI 
studies of image size evolution. It is not clear from the presented dataset whether we have seen emission from the counter jet or not. 
We predict that the Low Frequency Array will be able to observe the afterglow of GRB 030329 and many other radio afterglows, 
constraining the physics of the blast wave during its non-relativistic phase even further. 

Key words, gamma rays: bursts - radio continuum: general - radiation mechanisms: non-thermal 



^ 1. Introduction this burst was 5.5 xlO^^ergcm^^ in the 7-30 keV band, and l.lx 
H 10 '*ergcm"^ in the 30-400 keV band. The burst was followed 
Ci GRB 030329 has been a very distinctive event in many respects. ^ ^„ extremely bright X-ray afterglow, 1 .4 x 10-"'erg cm-^ s"' 
Residing at a redshift of 0.1685 dGreiner et al.| | 2003|) , i.e. at a jn the 2-10 keV ban d, detected by RXTE ~ 5 h after the burst 
luminosity distance of 802 Mpc (adopting a flat universe with ([Marshall & SwanklllOOl . The optical afterglow wa s detected, 
Qm = 0.27, Qa = 0.73 and Ho = 71kms Mpc ), it is one ^7 ^^^^^^^ ^^^^^ ^^^^^^ ^ ^^^^ ^2.4 mag (ISato et al.l 
atZT. ° ^y^^of ^".^fterglow has been observed 2003). A bright radio afterglow of 3.5 mJy at 8.46 GHz was de- 
(GRB 980425 at z = 0.0085 remains the nearest of the GRBs). tecte d by die Very Lar ge AiTay (VLA) on 2003 March 30.06 
GRB 030329 displayed one of the brightest afterglows ever, en- dBerger et al.ll200l . Around 7 days after the burst the op- 
abling the study of its evolution for a long time and in detail over ^-^^^ ^^^^^^^ ^1^^^^^ the signatur e of u nderlying sup ernova 
a broad range of frequencies, from X-ray to centimetre wave- sN 2003dh (iHiorth et all I200I IStanek et al.Fl200l . 
lengths. The afterglow of this GRB is one with the longest fol- ^j^^ afterglow was subsequently foUowed at X-ray, optical, mil- 
low ups ever, still visible in radio waves three years after the ^.^^-^ wavelengths, providing the exti-emely rich 
burst trigger It was also the first GRB to have provided unam- ^^^^^.^^ coverage of the tr ansient in a U die w avebands (e.g. 
biguous evidence of the long suspected (e.g. Galama et a . 1998) ^-^^^^ 200 4 HBTe t al. 200 ^ ^3^a bel et al. Mt 



association betwe en GRBs and supernovae (^ Hjorth et al. ,, 2003 . ; ighefli et al.1 12003t iKuno etaLi |2004; Ber geretaLi u2003). The 

I Stanek et alj | 2003[) . proximity of the GRB and die exceptional radio brightness of its 

GRB 030329 was detected and localised by the HETE-II afterglow made it possible to resolve the afterglow image with 

satellite ( Vanderspek et a l. 2003) on 29th March 2003, at UT ^ Very Long Basel ine Interferometi-y ( VLBI) campaign from 

11:37:14.7 and lasted mor e than 100 s. The measured fluence for g ^^^^ jTavlor et al.ll2004 up to 806 days (iPihlstrom et al.ll2007h 

Send offprint requests to: A.J. van der Horst, e-mail: ^^^^^ burst. 
avdhorst@science . uva . nl 
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Extensive radio follow up of the a fterglow h a s bee n 



reported ear lier by (iBerger et al.' '20031), 



• py ^ 

[Van der Horst et alj (l2005l) and 'Resmi et al 



Frail et all (l200l . 
(|2005|). These re- 



ports cover a period of up to 1 year post burst. In this paper, 
we report further extension of the low frequency radio follow up 
of the GRB 030329 afterglow, up to 1128 days after the burst, 
using the Westerbork Synthesis Radio Telescope (WSRT) and 
the Giant Metrewave Radio Telescope (GMRT). GRB 030329 is 
the first afterglow to be detected below an observing frequency 
of 1 GHz, and even as low as 610 MHz. In Section|2]we describe 
the observations and data analysis. In Section [3] we present the 
results from modeling the light curves of all the available radio 
data of the GRB 030329 afterglow, focusing on the blast wave 
physics in the non-relativistic phase of its evolution. In Section|4] 
we compare our modeling results with previous light curve stud- 
ies and with the VLBI image size evolution. Furthermore, we 
put constraints on the emission of the counter jet, and we show 
model predictions for radio afterglow observations with the Low 
Frequency Array (LOFAR). Section|5]summarises our results. 

2. Observations & Analysis 

GRB 030329 was observed with the WSRT and GMRT from 
325 MHz to 8.46 GHz, spanning a time range of 268-1128 days 
after the burst. The afterglow was clearly detected at all frequen- 
cies except for 325 MHz, where we only obtained upper limits. 
GRB 030329 is the first afterglow to be detected at frequencies 
less than 1 GHz: at 840 GHz with WSRT and even as low as 
610 MHz with GMRT. Here we describe the data reduction and 
analysis of our observations; the observational results are sum- 
marised in Table [U and Table ID for WSRT and GMRT respec- 
tively. 

2.1. WSRT Observations 

The first observations of GRB 030329 with WSRT were car- 
ried out at ~ 1.4 days after the burst, at 1.4, 2.3 and 4.8 GHz. 
After the clear detections of the afterglow, we started an inten- 
sive monitoring campaign at these three frequencies, o f whic h 
the first results were presented in IVan der Horst et al] (l2005l) . 
Here we present the results of follow up measurements up to 
1128 days after the burst in a wider frequency range, adding 
detections at 840 MHz and 8.4 GHz, and upper hmits at 3 50 
MHz. We used the Multi Frequency Front Ends ( lTanlll99ll) in 
combination with the IVC-hDZB backencQ in continuum mode, 
with a bandwidth of 8x20 MHz. Gain and phase calibrations 
were performed with the calibrator 3C286, though sometimes 
3C147 or 3C48 were used. The observations were analysed 
using the Multichannel Image Reconstruction Image Analysis 
and Display (MIRIAD) software package; except for the 350 
MHz observations, which were analysed using the Astronomical 
Image Processing System (AIPS). Table [1] lists the log of the 
observations, spanning a time range of 268-1128 days after the 
burst. We checked our results for consistency by measuring the 
flux of several nearby point sources in the primary beam, which 
were assumed to be constant in time. We found no significant 
correlated flux-variations between the various epochs. 

2.2. GMRT Observations 

The radio afterglow of GRB 030329 was first detected by GMRT 
at 1280 MHz on the 31st of March 2003, 2.3 days after the burst 



See Section 5.2 at |http://www.as tron.nl/wsrt/wsrtGuide/node6.htmll 



(iRao et al.ll2003h with a flux of 0.25 mJy. The afterglow was ob- 
served since then at 1280 MHz, 610 MHz and 325 MHz. We 
observed the afterglow at a total of 27 epochs (12 epochs at 
1280 MHz, 13 epochs at 610 MHz and 2 epochs a t 325 MHz) 
excluding the first year observations reported in iResmi et"al] 
(l2005l) (9 epochs at 1280 MHz). We have used a bandwidfli of 
32 MHz for all these observations. One of the three possible flux 
calibrators, 3C48, 3C147 or 3C286, was observed at the begin- 
ning and end of each observing session for about 15 minutes, as 
a primary flux calibrator to which the flux scale was set. Radio 
sources 1 125-H261 and 1021+219 were used as phase calibrators 
at 1280 MHz and 610 MHz, respectively. The phase calibrator 
was observed for about 6 minutes before and after an observing 
scan of about 30 to 45 minutes on GRB 030329. The data thus 
recorded were then converted to FITS format and analysed using 
AIPS. Fluxes of the individual sources were measured using the 
task 'jmfit' in AIPS. 

During the data reduction, we had to discard a fraction of the 
data affected by radio frequency interference and system mal- 
functions. At the end of this we could use an effective band- 
width of about 13 MHz, on average, for all the observations 
reported here. Estimated fluxes of other sources in the field of 
GRB 030329, i.e. in the primary beam, showed a variation of up 
to -20 % between observations at different epochs. A part of this 
variation appears to be correlated between different sources. This 
indicates the presence of possible systematic uncertainties in the 
flux calibration arising from instrumental effects, for which we 
attempt to correct in the following manner At each frequency 
band, we have chosen a reference observation: at 610 MHz we 
have chosen the field observed on September 2 2005 as a refer- 
ence, and at 1280 MHz the one observed on July 1 2005. For ev- 
ery other observation, we selected four sources in the field within 
5' of the GRB 030329 position and estimated their fluxes. Ratios 
of these estimated fluxes to those at the reference epoch were 
computed and averaged for each observation. The flux of the af- 
terglow measured at each epoch was then scaled using this aver- 
age flux ratio. The correction to the GRB flux obtained from this 
procedure was quite small, typically less than the map r.m.s. at 
most epochs, except at two instances at 1280 MHz where the cor- 
rection amounted to ~ 2cr. This gives us confidence that the esti- 
mated flux of the GRB is quite secure at both 610 and 1280 MHz. 
The final fluxes, with the corrections mentioned above incorpo- 
rated, are presented in Table |2] 

3. IVIodeling Results 

We have modeled the light curves o btained together with pre - 
viously report ed fluxes fro m WSRT ( Van der Horst et aLll2005h . 
GMRT (Resm i et al .ll2005h. and VLA, ATCA & Ryle Teles cope 
(IBerger et al.| [2003; Fra il et"ani2005t iPihlstrom et al.ll2007h . see 
Figure [T] The light curves show the characteristics that are ex- 
pected for the low-frequency part of a GRB afterglow. Since both 
the peak frequency v,n of the spectral energy distribution and the 
synchrotron self-absorption frequency Va are situated above the 
radio regime at early times and they both decrease in time, the 
light curves rise. When the peak of the broadband synchrotron 
spectrum, either or Va, has moved through the observing band, 

the light curves turn over and decline st eeply. 

According to the blast wave model (iRees & Meszarosll 19921; 
iMeszaros & Reeslll997l; IWiiers et al.1ll997l) the afterglow of a 
GRB is due to non-thermal synchrotron radiation emitted by 
shock accelerated electrons. The large amount of energy released 
during the burst, with a collimated outflow to start with, drives 
a powerful relativistic blast wave. As the blast wave propagates 
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Observing Dates 




Integration time 


Frequency 


Flux* 




(days) 


(hours) 


(GHz) 




2003 Dec 22.926 - 23.196 


268.577 


6.5 


2.3 


1618 + 31 


2003 Dec 25.005 - 25.188 


270.612 


4.4 


1.4 


2502 ± 139 


2004 Jan 25.833 - 26.333 


302.599 


12.0 


0.35 


< 951(3o-:i;. 


2004 Jan 29.822 - 29.961 


306.408 


3.3 


1.4 


1552+ 111 


2004 Jan 29.995 - 30.134 


306.580 


3.3 


4.8 


< 648 Oo--;^ 


2004 Jan 30.168 - 30.307 


306.753 


3.3 


2.3 


1389 ± 67 


2004 Jan 31.012 -31.137 


307.591 


3.0 


8.4 


815 + 225 


2004 Feb 10.790- 11.289 


318.555 


12.0 


0.84 


2332 ± 288 


2004 Mar 26.667 - 26.934 


363.816 


6.4 


4.8 


597 ± 27 


2004 Mar 27.664 -28.152 


364.424 


11.7 


1.4 


1318 ± 104 


2004 Mar 28.862 -29.153 


365.524 


7.0 


2.3 


871 ± 29 


2004 Apr 11.623 - 11.880 


379.267 


6.2 


0.84 


1525 + 389 


2004 May 2.566 - 3.065 


400.331 


12.0 


0.35 


< 1305 {3cry 


2004 May 19.634-20.018 


417.342 


9.2 


1.4 


1824 + 100 


2004May 22.511 -22.794 


420.168 


6.8 


2.3 


933 + 34 


2004 Jul 3.396 - 3.729 


462.078 


8.0 


2.3 


707 + 39 


2004 Jul 4.394 - 4.726 


463.076 


8.0 


4.8 


329 + 27 


2004 Aug 1.317- 1.598 


490.974 


6.7 


1.4 


622 + 95 


2004 Sep 25.328 - 25.667 


546.013 


8.1 


4.8 


274 + 34 


2004 Nov 2.063 - 2.271 


583.683 


5.0 


2.3 


543 + 46 


2004 Nov 11.039 - 11.371 


592.721 


8.0 


1.4 


1162 + 63 


2004 Nov 12.036 - 12.536 


593.802 


12.0 


0.84 


1306 + 366 


2004 Nov 20.014 -20.486 


601.766 


11.3 


0.35 


< 2124 {3cry 


2005 Mar 24.673 -25.172 


726.439 


12.0 


0.84 


1199 + 238 


2005 Mar 26.667 -27.167 


728.433 


12.0 


0.35 


< 996 (3fr)'- 


2005 Apr 9.629- 10.129 


742.395 


12.0 


1.4 


1078 + 72 


2005 Apr 10.627 - 11.126 


743.892 


12.0 


2.3 


504 + 48 


2005 May 14.534 - 14.919 


777.242 


9.2 


4.8 


409 + 27 


2005 May 15.531 - 15.883 


778.223 


8.1 


8.4 


< 309(3o-:i^ 


2005 Nov 27.993 - 28.493 


974.759 


12.0 


4.8 


150+19 


2005 Dec 7.004 - 8.292 


984.163 


13.8 


2.3 


318 + 58 


2005 Dec 9.000 - 10.233 


986.132 


12.5 


1.4 


842 + 127 


2006 Apr 9.637 - 10.129 


1107.399 


11.8 


0.84 


817 + 392 


2006 Apr 30.573 - May 1.072 


1128.338 


12.0 


4.8 


157 + 22 



" The indicated time is the logarithmic average of the start and end of the integration. 
* The measurement uncertainties are given at Icr level. 

A meaningful formal flux measurement at the GRB position cannot be determined because of confusion by a very nearby, bright source. 

Formal flux measurement at the GRB position gives 1+216 yujy. 

Formal flux measurement at the GRB position gives -69 + 103 fiJy. 



into the circumburst medium, the electrons in the medium are 
accelerated to a power-law distribution of Lorentz factors with 
power-law index p. These relativistic electrons gyrate in the 
post-shock magnetic field and emit synchrotron radiation which 
is seen as the afterglow of the GRB. The power-law distribu- 
tion of the electrons results in a power-law spectrum of the af- 
terglow. Meanwhile the shock wave decelerates as it propagates 
into the circumburst medium. Assuming mass-energy conserva- 
tion across the shock front, it has been shown that the Lorentz 
factor of the blast wave falls oflF as a power-law with its radius. 

It is expected that the decelerating shock wave becomes non- 
relativistic after a few weeks. Somewhat earlier, the sideways 
expansion of the initially tightly collimated outflow would be- 
come important, and by the time of the non-relativistic transition 
the shock wave becomes nearly spherical. In the non-relativistic 
phase the evolution of the shock wave can be described using 
the Sedov-Taylor self sim ilar solutions. A d etailed description of 
this phase can be found in'F rail et all ( 12000') (hereafter FWKOO). 
Observations of the broadband afterglow during this phase can 
be used to estimate various physical parameters related to the 
explosion, such as the amount of energy released during the ex- 



plosion, fractional amount of energy in the accelerated electrons 
and in the post-shock magnetic field, and the density and struc- 
ture of the circumburst medium. These parameters can also be 
estimated, independently, by modeling the evolution of the after- 
glow in the relativistic phase, but this suffers from uncertainties 
related to collimation geometry and relative orientation of the 
observer, problems that do not plague the non-relativistic phase. 

The broadband afterglow of GRB 030329, from radio to X- 
ray frequencies, can be modeled by the standard rel ativistic blast 
wave model, assuming either a ' double jet model' dBerger et al.l 
l2003i) or a 'refreshed jet model' (iResmi et alj|2005l) . In both flie 
models, the early-time optical and X-ray light curves are ex- 
plained by a jet with a small opening angle of ~ 5°. The double 
jet model assumes a co-aligned wider jet component (~ 20°) that 
caiTies the bulk of the energy and produces the later time light 
curves. In the refreshed jet model, the initial jet is re-energised 
by the central engine during its lateral expansion, re-collimating 
it to a wider opening angle. This refreshed jet then produces 
the late-time emission. Unfortunately, we cannot distinguish be- 
tween the two models with all the available data, including the 
data presented in this paper. 



4 



Van der Horst et al.: GRB 030329 radio afterglow deep into the non-relativistic phase 



Table 2. Log of GMRT observations of GRB 030329. 



Observing Dates 


AT 


Frequency 


Flux* 




(days) 


(MHz) 


(mJy" 


2004 Jan 4.97 


281.49 


610 


0.90 ± 0.20 


2004 Feb 26.63 


334.15 


610 


1.14 ±0.25 


2004 Apr 10.49 


378.01 


610 


0.61 ±0.11 


2004 Jun 27.28 


455.80 


610 


1.01 ±0.18 


2004 Jun 28.28 


456.80 


1280 


1.11 ±0.13 


2004 Jul 18.31 


476.83 


1280 


1.04 ±0.11 


2004 Jul 20.28 


478.80 


610 


0.84 ±0.17 


2004 Aug 20.04 


509.56 


1280 


1.03 ±0.11 


2004 Sep 06.00 


526.52 


610 


0.98 ±0.15 


2004 Sep 17.96 


538.48 


1280 


0.48 ±0.10 


2004 Dec 3.00 


614.52 


610 


1.08 ±0.35 


2005 Jan 6.82 


649.34 


610 


0.50 ±0.11 


2005 Feb 8.81 


682.33 


325 


< 1.5(3o-:i^ 


2005 Feb 13.84 


687.34 


1280 


0.97 ±0.14 


2005 Feb 21.71 


695.23 


610 


0.47 ±0.13 


2005 Mar 10.77 


712.29 


1280 


0.72 ± 0.07 


2005 Mar 18.66 


720.18 


610 


0.69 ±0.13 


2005 Jun 17.54 


811.06 


610 


< l.l(3cr:j^ 

< 1.8(30-)'' 


2005 Jun 28.59 


702.11 


325 


2005 Jul 1.19 


824.71 


1280 


0.53 ± 0.09 


2005 Sep 2.99 


888.51 


610 


0.68 ±0.13 


2005 Oct 8.20 


923.72 


1280 


0.49 ± 0.09 



* The measurement uncertainties are given at lo" level. 

* A meaningful formal flux measurement at the GRB position cannot be determined because of confusion by a very nearby, bright source. 
Formal flux measurement at the GRB position gives 0.74 ± 0.36 mJy. 



After 80 days the observed radio light curves flatten, which 
can be explained by a transition into the n on-relativistic phase 
of the blast wav e jVan der Horst et aL 2005; Resmi et al. 2005; 
iFrail et al.ll2005h . It was suggested in lVan der Horst e t al. (2005) 
that this late-time behaviour could also be explained by a third 
jet-component with an even wider opening angle than the first 
two. The latter model, however, is excluded by the observa- 
tions below 1 GHz (see Figure [U, which leaves the model 
with the non-relativistic phase after 80 days as the preferred 
model for the late-time behaviour of the blast wave. The precise 
value of the time of the transition into the non-relativistic phase 
varies between 60 to 80 days because of the different model- 
ing methods applied by various authors. The transition into the 
non-relativis t ic pha se is also observed at X-ray frequencies by 
iTien^o et al.l (12004 *). who estimate a transition time of ~ 44 days. 
The X-ray light curve, however, is less well sampled than the 
radio light curves at centimetre wavelengths, which makes the 
determination of the transition time uncertain. 

Because most of the GRBs occur at cosmological distances, 
they are not bright enough to be observed at late times, and in 
fact most of them fade below the detection limits much before 
the start of the non-relativistic phase. GRB 030329, being one 
of the closest GRBs, provides us with a unique example of late 
time observations deep into the non-relativistic phase. Since the 
first months of the radio afterglow have been extensively dis- 
cussed in the literature, we will focus on the late-time behaviour 
of the light curves. Here we present a detailed analysis of the 
non-relativistic phase of the afterglow. 



3.1. Spectral & Temporal Slopes 

We investigate the spectral and temporal slopes in the non- 
relativistic phase by performing a joint temporal and spectral 
power-law fit of F,, = F„ • (vghz/4.86)-'' ■ (fdays/lOO)"". Because 



of the uncertainty in the exact transition time into the non- 
relativistic phase, fNR, and because the value for the temporal 
slope could be altered by the smoothness of the transition, we 
use all the available data after 100 days at 2.3, 4.8, 8.4, 15 and 
22.5 GHz. At 1.28 & 1.4 GHz we use all the data after 300 days, 
as the light curve at these frequencies peaks around 100 days. We 
find a = 1 .08 + 0.03, ;S = 0.54 + 0.02 and F„ = 2.69 + 0.06 mJy. 

After 100 days the peak frequency and synchrotron self- 
absorption frequency have passed through the observing bands 
we used for this fit, and we do not find any evidence for a 
chromatic break, caused by the so-called cooUng break, up to 
1128 days. So in this case the spectral slope /? is given by 
ip - l)/2, and the temporal slope a is given by 3{5p - 7)/10 
or (7p - 5)/6, for a homogeneous and wind circumburst medium 
respectively. The power-law index p of the electron energy dis- 
tribution is thus given by 2/3+1, and (10a-i-21)/15 (homogeneous 
medium) or (6a -i- 5)/7 (wind medium). For our fitted values we 
find p = 2.09 + 0.03 from the spectral slope, and p = 2.12 + 0.02 
(homogeneous) or p = 1.64 + 0.02 (wind) from the temporal 
slope. This indicates that the circumburst medium is homoge- 
neous in the non-relativistic phase of the blast wave evolution, 
and that a wind medium is completely rejected in this phase. If 
one assumes that the circumburst medium density is a power-law 
of the radius r with index k, i.e. proportional to r"*, the value of 
k and its uncertainty can actually be calculated directly from a 
andjS: k = (5a'-15j6-t-3)/(a'-4j6-H2).Forourfitted values we get 
k - 0.33;^Q4|, consistent with a homogeneous medium (k = 0); 
and this value is inconsistent with a wind medium (k - 2) at the 
5cr level. 

The value of /:> = 2.1 is in agreement with values found by 
other authors, although it does matter which data one includes, 
i.e. whether one takes only the data after 100 days or also earlier 
data. Including earlier data gives higher values for a and /3, and 
as a result a higher value for p, with a significantly higher re- 
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Fig. 1. Modeling results of the afterglow of GRB 030329 at centimetre wavelengths. Our light cur ves obtain ed with WSRT an d 
GMRT are shown together with previously reported fluxes from WSRT ( Van der Horst et all 12005 !). GMRT dResmi et al.ll2005h . 
and VLA, ATCA & Ryle Telescope (Bergeret al. 2003; Frail et al. 2005). The filled circles and squares are WSRT and GMRT 
measurements respectively; the open diamonds are VLA, ATCA & Ryle Telescope measurements; the open triangles are 3cr upper 
limits. Three fits to the data are shown: the dotted line represents a fit to the first 100 days of radio observations with a wide jet 
expanding in a homogeneous medium; the solid line corresponds to a model in which the blast wave becomes non-relativistic after 
80 days; the dash-dotted line corresponds to a model in which a third jet-component with an even wider opening angle is present. 
The latter model is excluded by the observations below 1 GHz, which leaves the model with the non-relativistic phase after 80 days 
as the preferred model for the late-time behaviour of the blast wave. 



duced chi-squared. This indicates that if one includes part of the 
data before 100 days, a steeper decay phase is also sampled, i.e. 
the jet-spreading phase. And thus our choice of fitting the data 
after 100 days is a valid assumption. 



3.2. Physical Parameters 

We have used 1128 days of broadband radio observations 
(610 MHz to 43.3 GHz) of the afterglow of GRB 030329 to 
model the dynamical evolution of the afterglow as well as to 
constrain the explosion energy. In Figure [T] we compare our 
m odel predictions wit h the observations. The model presented 
m IWn der Horst et al ] (l2005i) . which fitted the data up to eight 
months, per fectly fits the data up to m ore than three years after 
the burst. In I Van der Horst et al.l ( 12005 ) values for the peak fre- 
quency Vn, of 35 GHz and for the self-absorption frequency 
of 13 GHz were found, both measured at the jet-break time of 
10 days; the flux at v,n at that time was 61 mJy. After about 80 



days the afterglow shows a flattening of the light curve, which is 
the start of the non-relativistic phase of the explosion: fNR- 

The fitted spectral parameters may now be used to derive the 
physical parameters of the explosion. The break frequencies and 
the peak flux, estimated deep in the non-relativistic phase, i.e. 
at a reference time to » Tnr, can be used to yield an estimate of 
the blast wave energy Esj and the ambient baryon density n,, us- 
ing the Sedov-Taylor solution for the blast wave, in the manner 
adopted by FWKOO for GRB970508. Two other physical param- 
eters that determine the evolution of the radiation are the fraction 
of total energy in relativistic electrons (se) and in the post-shock 
magnetic field (sb). In order to determine these four quantities, 
one requires the measurement of four spectral parameters, tradi- 
tionally the three break frequencies and the flux normalisation. 
In the late phase, however, direct determination of the cooling 
frequency is difficult, since the afterglow is not detectable at fre- 
quencies above radio bands. We therefore express the physical 
parameters as a function of the ratio e, = Se/sB; £r = 1 would 
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signify an equipartition of energy between the magnetic field and 
the relativistic particles. 

Following Eq. (5) of FWKOO, we may then write the post- 
shock electron number density as 



Bl 



(1) 



We insert this in Eq. (A 12) of FWKOO, and invert their 
Eqs. (AlO) to (A12) while making use of the relations (A6) to 
(A8). This yields the following set of expressions for the FWKOO 
model parameters: 
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In the above, ro is the radius of the blast wave, Bq is the post- 
shock magnetic field, and yo is the lower cutoff of the power-law 
distribution of relativistic electron Lorentz factors, all at the ref- 
erence time fo- From these, the blast wave energy may be com- 
puted as: 



Est = noup 
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In Eq. ^ through (|7]i, Vao.GHz and Vmo.GHz are the two break fre- 
quencies, in GHz units, at the reference time fo - fo.d days. The 
expressions assume that Vao.GHz > VmO.GHz, which is indeed the 
case in GRB030329 at these late times. Fmo.mjy is the normalisa- 
tion of the flux at an observing frequency v » Va, Vm, expressed 
as 



Fy{t = fo) mjy = F,nO,mJy 



-(p-l)/2 



(8) 



z is the redshift of the burst and ^28 the corresponding luminosity 
distance. The thickness of the post-shock emitting region at any 
time is assumed to be r/?/, where r is the radius of the blast wave 
and 77 = IOt/io. fi and /3 are integrals over the synchrotron func- 
tion defined in FWKOO; both are functions of p. The quantity 



^, close to unity, is an equation of state-dependent normalisa- 
tion factor for the blast wave radius (FWKOO). represents the 
mass fraction of hydrogen in the circumburst medium, and may 
nominally taken to be 0.75. 

Evaluating spectral parameters from the fitted model, one 
finds, at a reference time fo - 500 days: VmO.GHz = 2.56 x lO""*, 
Vao.GHz = 0.498 and Fmo.mjy = 28.4. Using z = 0.1685 for 
GRB030329 and the fitted value of p = 2.1, one then estimates 



Est = 2.6 x 10^' erg 



Hi = 0.35 cm"^ (-^] 
\'7io/ 



£1- 

0.41 



0.12 



(9) 
(10) 



The blast wave radius at 500 days works out to be ro = 
0.49(ei/77io)"''''^ pc. The corresponding post-shock magnetic 
field is Bo = 0.036(ei//7io) G, and the lower cutoff of 
the electron Lorentz factor distribution at that time is yo = 
5.5(e,/77io)° '2. These yield fie = firfiB = 0.085(er/?7io)"-2'*. 



4. Discussion 

The data presented in this paper give an unprecedented view of 
the non-relativistic evolution phase of a GRB blast wave, be- 
cause of the wide range covered in both frequency and time. This 
gives us the opportunity to compare the physical parameters that 
we have derived from the very late-time data with the physi- 
cal parameters derived from the early-time data, when the blast 
wave was still extremely relativistic. From the emerging physi- 
cal picture we put constraints on the emission from the counter 
jet, and we compare our findings with the results from VLBI 
measurements of the source size evolution. 



4.1. Relativistic versus Non-Relativistic 



In Section 13.21 we calculated the total energy in the blast wave 
Est and the density of the circumburst medium «, as functions 
of the ratio = Se/sB, assuming that the blast wave was in 
its non-relativistic phase. These parameters are determined in 
IVan der Horst et alT(l2005h in the relativistic phase as functions 
of the cooling frequency, adopting Vc = 10'^ Hz. We calculate 
the energy and density in the relativistic phase as functions of 
fir, a better way to compare the parameters in the two different 
phases. Since the values of Se = 0.25vyj*3 and sb = 0.49 v^^j^'* in 
der Horst et al. I (l2005h indicate near-equipartition, the val- 
ues for the density and energy do not differ much from the values 
obtained there: Hr = 0.78xey^cm"^ £^450 = 2.6x 10^' xe^'' erg 
and Ei-^cor = 3.4 x 10^" x s^J* erg, the latter being the beaming- 
corrected energy. The cooling frequency can be found to be 
Vc - 1.6 X lO'"* X Ey^^ Hz, validating the assumptions in pre- 
vious studies on being ~ lO'-' Hz at 10 days. 

We find that the density derived in the non-relativistic phase 
is a factor of two smaller than the density derived in the rela- 
tivistic phase. The beaming-corrected energy iSr.cor, however, is 
a factor of ~ 7 times smaller than the total energy Est in the 
blast wave derived from the non-relativistic evolution. iFrail et al.l 
(i2005i) modeled the first year of observations with the VLA and 
ATCA, and found a total kinetic energy of 9.0 x 10^" ergs; 
using only the data before 64 days they derive an energy of 
6.7 X 10^*^ ergs (see also Berger et al. 2003), a factor of ~ 2 larger 
than our Ey coi-\ and using only the data after 50 days, to obtain 
an estimate for Est, they find an energy of 7.8 x 10^" ergs, a 
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factor of ~ 4 smaller than our value for Esj. In these models the 
range for the circumburst densit y is ~ 1 - 3 cm"-', a bit larger 
than our values. Granot et al. (2005|) also determine the energy 
and density with their models, and find a similar value for the 
collimation corrected energy as we do. Their value for the den- 
sity is, however, an order of magnitude larger, but they note that 
this can be attributed to the fact that the density depends strongly 
on the precise value of Va and v^. 

Given the differences in the methods used by different au- 
thors and the uncertainties in the assumptions made to estimate 
these numbers, they can all be considered to be comparable. It is 
not possible to make any definite statements about significant 
differences in energies derived from the relativistic and non- 
relativistic phase. If, however, the somewhat larger estimate of 
the total energy in the non-relativistic phase is indeed true, then 
two possible explanations may be advanced for this: either the 
beaming correction in the relativistic phase is too strong, giving 
a too small value for the beaming corrected energy; or Esj is 
over-estimated because of non-isotropy in the emission coming 
from the blast wave in the non-relativistic phase. In the latter 
case it could be that the blast wave is not completely spherical 
yet, but still the evolution is well described by the Sedov-Taylor 
solution; or that the blast wave is spherical, but the emission is 
not coming from the blast wave isotropically. In both situations 
the value of Est that we derived would be an over-estimate of 
the true value. 

The energies that we and iFrail et al.l (l2005i) derive indicate 
that our estimate of the relativistic beaming and the assumption 
that the blast wave becomes spherical at f ^ Jnr are valid. The 
latter is important for testing the models that describe the lat- 
eral spreading of the collimated outflow after the jet-break time, 
when the Lorentz factor drops below the inverse of the half- 
opening angle of the jet. Some (semi-analytical) models (e.g. 
[Rhoads 1999) assume a very rapid sideways expansion of the jet 
with a lateral expansion velocity of the order of the velocity of 
light, resulting in an exponential growth of the jet half-opening 
angle with radius. Hydrodynamical simulations, however, show 
a very modest degree of lateral expansion as long as the jet is 
relativistic (for an extensive review, see Granot 2007). In the lat- 
ter model the outflow is still strongly collimated when the blast 
wave becomes non-relativistic, while the first model predicts that 
the blast wave is (almost) spherical at f ^ Jnr, which is favoured 
by our analysis, since there is no significant change in temporal 
slopes well after Jnr- 

4.2. Counter Jet Emission 

It is expected that there are two collimated outflows formed at 
the collapse of a massive star: two jets pointed away from each 
other. When the counter jet becomes non-relati vistic, the emis- 
sion is no longer strongly beamed away from us. lGranot & Loebl 
(l2003h predicted a re-brightening in the radio light curves when 
that occurs, altho ugh it would be difficult to detect in the case 
of GRB 030329. Li & Song! (|2004 calculated that the observed 
time at which the counter jet becomes non-relativistic is 5 times 
larger than Jnr, because of light travel time effects. They claim 
that the observed flux coming from the counter jet at that time 
is comparable to the flux of the jet coming towards us at Jnr. 
This results in a rapid increase in flux with a peak at 5 x fNR, 
and after the peak the light curve declines steeper than before 
the rapid increase until it relaxes to the original light curve be- 
haviour. Li & Song (2004) predicted such a feature in the ra- 
dio light curves of GRB 030329 at ~ 1.7 years after the burst. 
From Figure[T] however, it is clear that this is not observed up to 



3 years after the burst (which was also noted, with ob servations 
up to 2 years after the burst, bv'P ihlstrom et alj|2007h . 

The apparent discrepancy between the predictions for the 
counter jet and the observed light curves can be explained 
by looking carefu lly at the assumptions. The calculations by 
iLi & Sona (|2004|) are valid for those observing frequencies at 
which the light curve peaks when the blast wave is still ultra- 
relativistic and narrowly collimated. This means that it can only 
be applied to the light curves at 8.4 GHz and higher frequencies. 
Adopting our value for Jnr of 60-80 days, this would mean that 
the 8.4 GHz light curve has to rapidly increase to a peak of 1- 
2 mJy at 300-400 days (5 x Jnr). This is not observed, although 
it is hard to make definite statements about this, since there are 
only observations up to 360 days. 

The light curves at 2.3 GHz and lower frequencies are not 
expected to show this kind of re-brightening, because they peak 
at or after Tnr- This results in a flattening of the peak to a width 
of 2 X fNR. Looking at the low frequency light curves this could 
be the case, although the scatter in the data is too large to make 
any quantitative statements about this. 

The best constraints on the c ounter jet could co me from the 
4.8 GHz light curve. Following iLi & Sond ( |2004 the flux at 
300-400 days has to increase to 4-6 mJy, while the measured 
flux at that time is 0.6 - 0.8 mJy. This means that the observed 
flux coming from the counter jet is suppressed by a factor of 
~ 8, while at 8.4 GHz the limit was a factor of ~ 4. The flux 
from the counter jet, however, is probably lower than that pre- 
dicted bv lLi & Son3 (|2004 . because, especially at 4.8 GHz, the 
light curve peaks later than the jet-break time. This means that 
the two jets, which already have quite a large opening angle of 
~ 20° to start with, cannot be treated as two narrowly collimated 
outflows anymore because of lateral spreading of the jet. 

Concluding, we cannot say, from the light curves presented 
here, whether we have seen emission from the counter jet or not. 
The fact that we do not see a late-time re-brightening at high 
radio frequencies could be due to the fact that the outflow is not 
very narrowly collimated. A flattening of the peak of the light 
curves at low radio frequencies, caused by the emission coming 
from the counter jet, could be present, but the scatter in the data 
prevents us from drawing any firm conclusion regarding this. 

4.3. Source Size Evoiution 

VLBI observations of GRB 030329 make it possible to study 
the image size of the afterglow. The full set of VLBI mea- 
suremen ts up to 806 days after the burst h as been presented in 
IPihlstron TetaLl (|2007|) . iGranot et all jlOOS ) have discussed the 
implications of the measured image sizes up to 83 days after 
the burst, which have mainly remained the same after includ- 
ing the latest measurements. They conclude that a homogeneous 
medium gives a better fit to the image size evolution than a stel- 
lar wind environment, although it is hard to rule out the latter 
due to the measurement uncertainties. Furthermore, they try to 
constrain the amount of lateral spreading, i.e. whether the lateral 
expansion velocity is of the order of the velocity of light, or that 
there is hardly any lateral spreading until the blast wave becomes 
sub- or non-relativistic. This is, however, hard to constrain, again 
because of the uncertainties in these source size measurements. 

A n important paramet er th at comes out of t he VL BI studies 
is fNR. lGranot et al.l(l2005h and lPihlstrom et al.l(l2007l) claim that 
fNR ~ 1 year from looking at the evolution of the apparent ex- 
pansion velocity, and from theoretical model fits to the source 
size evolution. The models in which there is rapid lateral expan- 
sion, their model 1, fNR ~ 1 year is indeed obtained, but the 
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jet-break time in that case is f ound to be ~ 1 month instead of 
10 days (e.g., see Figure 4 of' Pihlstrom et al.ll2007l) . Since the 
jet-break time is quite well determined from the light curves, the 
value of fNR in this model 1 is rather uncertain. The models in 
which there is no lateral expansion until the blast wave becomes 
non-relativistic, their model 2, the non-relativistic transition ac- 
tually happens at 60-80 days after the burst; the blast wave then 
becomes spherical on a time scale of 1-3 years. So the fits of 
model 2 do indeed indicate a value of Jnr ~ 60 - 80 days. 

From the evolution of the average apparent expansion ve- 
locity the value of Tnr can also be deduced. The uncertainties 
and correction factors in deriving the average velocity are dis- 
cussed thoroughly in Pi hlstrom et al. (2007 ). From their Figure 
3 it seems that the transition happens at ~ 100 - 200 days, when 
the average apparent expansion velocity is ~ 2 times the speed 
of light. The uncertainties in the correction factors applied, how- 
ever, make it hard to state that this is significantly higher than the 
value of fNR that is deduced from the light curves. The issue that 
still remains, is that the value of Jnr that various authors have 
estimated from energy considerations is a facto r of a few higher 
than 60-80 days (e.g. Granotetal. 2005; Va n der Horst et all 
12005 ). This is, however, a very rough estimate, and a discrep- 
ancy within a small factor should not be regarded as disturbing. 
In conclusion, the value of Jnr determined from the VLBI source 
size measurements is not at significant odds with the value of Jnr 
derived from light curve studies. 

4.4. Low Frequency Array 

From our analysis of the broadband radio afterglow of 
GRB 030329 we can calculate predicted light curves for future 
radio telescopes operating in the metre wavelength regime. As 
an example we explore the possibilities for the Low Frequency 
Array (LOFAR) to detect the GRB 030329 afterglow and other 
GRB afterglows. 

LOFAR will be a major new multi-element, interferomet- 
ric, imaging telescope designed for the 30-240 MHz frequency 
range. LOFAR will use an array of simple omni-directional an- 
tennas, whose electronic signals are digitised, transported to a 
central digital processor, and combined in software to emulate 
a conventional antenna. LOFAR will have unprecedented sensi- 
tivity and resolution at metre wavelengths, and will be the first 
of a new generation of radio telescopes to become fully oper- 
ational, i.e. in early 2009. This sensitivity and resolution will 
give the GRB community the opportunity to study bright after- 
glows on even longer timescales than with observations at cen- 
timetre wavelengths. For a concise description of LOFAR and 
the Transients Ke y Project, in wh ich GRBs are among the prime 
targets, see .Rottgering et al.l (12006) and Fender et al. (2006). 

The fact that our e arly-time model calculations from 
der Horst et al. I (120051) gave such good predictions for the 
late-time behaviour, in particular at observing frequencies below 
1 GHz, gives us confidence to extrapolate the modeling results 
of the radio afterglow of GRB 030329 to the LOFAR observ- 
ing range, see Figure |2] The predicted light curves show that 
GRB 030329 will be observable in the high band of LOFAR 
(120-240 MHz), but not in the low band (30-80 MHz). The Hght 
curves are expected to peak in 2009, when LOFAR will be fully 
operational, and even later going down from 240 to 120 MHz; 
the expected peak flux is ~ 0.3 - 0.4 mJy in the LOFAR high 
band. We also calculated light curves for GRB 030329 if it were 
situated at a redshift of 1 instead of 0.16. The resulting fainter 
afterglow can also be detected, although with longer integration 
times, i.e. on the order of a day instead of an hour 
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Fig. 2. The predicted light curves of GRB 030329 at three fre- 
quencies within the LOFAR observing range. The solid line 
and dash-dotted line correspond to the two models shown in 
Figure [Tl the dashed line shows the predicted light curve of 
GRB 030329 when situated at a redshift of 1 instead of 0.16. 
The horizontal dotted lines are the 1-cr sensitivity limits after 
4 hours of observing with a bandwidth of 4 MHz. It shows that 
it is quite feasible to detect GRB afterglows similar in brightness 
to that of GRB 030329 with LOFAR on timescales of months to 
decades; and that fainter afterglows can also be detected after 
long integration times. 



5. Conclusions 

We have presented a detailed study of the late-time radio af- 
terglow of GRB 030329. We have obtained measurements with 
the WSRT and GMRT, spanning a spectral range of 325 MHz- 
8.4 GHz and a temporal range of 268-1 128 days after the burst. 
Combined with all the already published radio observations of 
this afterglow, from WSRT, GMRT and other large radio tele- 
scopes, we have studied the physics of the blast wave in the non- 
relativistic phase and compared these results with those from 
studies of the relativistic phase. The well-sampled late-time light 
curves made it possible to determine the index of the electron en- 
ergy distribution accurately, and to confirm that the circumburst 
medium is homogeneous. 

The energy of the blast wave and density of the circumburst 
medium, determined from the non-relativistic evolution of the 
blast wave, are comparable to findings by several studies based 
on earlier time observations. We have shown that the blast wave 
is spherical, or almost spherical at least, at the time that the blast 
wave becomes non-relativistic: Tnr ~ 80 days. In contrast with 



Van der Horst et al.: GRB 030329 radio afterglow deep into the non-relativistic phase 



9 



some predictions, a radio re-brightening due to the counter jet 
becoming non-relativistic, is not observed. The existence of a 
counter jet cannot be ruled out, since it is possible that the peaks 
of the light curves at low radio frequencies are flattened due to 
this counter jet. We have also shown that the high-resolution 
VLBI measurements of the afterglow image size are in agree- 
ment with our light curve studies. In particular, the value of ?nr 
derived from modeling the image size evolution does not dif- 
fer from our findings significantly within the measurement and 
modeling uncertainties. 

The afterglow of GRB 030329 will be foUowed-up at ra- 
dio frequencies for many years to come. The current brightness 
and the fact that the flux drops logarithmically, make it possi- 
ble to continue studying this afterglow for the next decade at 
least. VLBI capabiUties will increase significantly in the coming 
years, so the image size evolution can also be observed. A new 
generation of radio telescopes will also be able to observe this 
afterglow. The first one of this new generation to come on-line 
is LOFAR, which will be able to detect GRB 030329 at metre 
wavelengths. Besides this extremely bright afterglow, LOFAR 
wiU also be able to study many afterglows at frequencies and on 
timescales that so far have been unexplored. 
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